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1 Introduction 

The inflationary scenario is now regarded as an indispensable ingredient in modern cosmology[l- 
5]. One of the most important predictions is the explanation of the origin of the primordial 
curvature perturbation as vacuum quantum fluctuations. The predicted curvature perturba- 
tion has definite statistical properties, such as the adiabaticity, nearly scale-invariance and 
nearly Gaussianity [6-10]. It turns out that these predictions are consistent with the current 
observation of cosmic microwave background (CMB) [11, 12]. Nevertheless, we still do not 
know the detailed model for the inflationary scenario. 

In this respect the possible observation of Non- Gaussianity in the primordial curvature 
perturbation in CMB temperature fluctuation has attracted much attention recently because 
it has a potentiality to discriminate various models of the inflationary scenario. To extract the 
useful information from the observational results it is necessary to have definite theoretical 
predictions in the nature of the Non-Gaussianity. This is our purpose in this paper. 

Here we shall focus on the so-called "local form" bispectrum [13-15] and trispectrum[16] 
in the primordial curvature perturbation. In general the bispectrum and trispectrum are 
defined as follows. 

3 

(QC(k*)) = (2^) 3 5 3 (^k,)S c (A: 1 ,A; 2 ,A;3), (1.1) 

i=l i 
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4 

<I]C(ki)> = (2vr) 3 5 3 (J]k i )r c (fc 1 ,fc 2 ,fc3,fc4). (1.2) 

i=l i 

The local form of them is given by the product of the power spectrum. The local form 
bispectrum and trispectrum are characterized by three parameters /nl, Ta/x an d 9nl as 
follows. 

B C = ^/nl [P c (fci)P c (fc2) + (2 perm.)] , (1.3) 

T c =r NL [Pc(|ki + k 3 |)P c (fc 3 )Pc(fc4) + (11 perm.)] 
54 

+ ^NL[Pc(fcl)Pc(^)Pc(fc3) + ( 3 P erm -)1- (1-4) 

There is an infinite number of higher moments defined as above. In this paper we focus only 
on f NL and r NL . The current best limit is f NL = 32 ± 21 (68% CL;[11, 12]). The Planck 
satellite is expected to reduce the error bar by a factor of four [15]. 

Recently, various consistency relations satisfied by non-linear parameters have emerged 
as a useful test of a large class of inflationary models. In the single field case, it is shown that 
a convincing detection of f^L ^> 1 or r ^> 1 rules out all single-field inflation models [17-25]. 
In the multi field case, Suyama and Yamaguchi proved the following inequality between f^L 
and tjsil [26] which we call SY inequality. 

TNL > ( g fNL J , (1-5) 

This inequality is derived using the 5N formalism [9, 27-30] in which the curvature pertur- 
bation is described by the departure of the e-folding and expanded by the scalar fields at 
near horizon crossing. In the original SY inequality, it is further assumed that non-linear 
parameters fjsiL and tnl are generated only on super horizon scales and the SN expansion 
is truncated up to the second order, that corresponds to the tree level in a diagrammatical 
method. As the extension of the SY inequality, Suyama et al. [31] pointed out that the 
truncated 6N expansion contains a part of 1-loop corrections and the original form of the 
SY inequality is still satisfied including such 1-loop corrections. In our previous work, we 
truncated the 5N expansion up to the forth order which corresponds to take into account all 
of the 1-loop corrections [32]. On the other hand, as pointed out in [33], the SY inequality 
is automatically satisfied by their definition of /nl and tnl regardless of the details of the 
inflationary model in infinitely large universe, and it is also suggested that the definition 
itself leads to the weaker condition for a finite observable universe 1 . Thus there remains 
still a possibility that the SY inequality is violated in a finite observable universe, and it is 
interesting to investigate the generality of the SY inequality from a physical point of view. 
This is important because future CMB observation such as the Planck satellite might detect 
these parameters in a good accuracy and give us a useful information on the inequality. 

x The Eq. (9) in [33] is given by 

>(§/*!,) (1.6) 

The second term of the Right Hand Side means the correction of a finite volume of the observable universe. 
This term vanish in the infinite universe, but always become negative in a finite observable universe. Therefore, 
this inequality is weaker than the SY inequality. However, it is not trivial what the specific difference between 
Eq. (1.6) and the SY inequality shown form the 5N formalism is 
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To study the generality of the SY type inequality, we need to extend the SY inequality 
up to all order of the SN expansion, that is all loop corrections because there is a possibility 
that higher order loop corrections give a large contribution in the evaluation of f^L and t^l- 
We address this question in this paper. We will show that the same expression with the 
original SY inequality is still satisfied approximately as far as the curvature perturbation 5N 
is expanded by linear scalar fields (namely, no vector fields). 

The organization of the present paper is as follows. In Sec. II, we explain some notations 
used in this paper. In Sec. Ill, we explain about the 5N formalism briefly. Then, we show 
the detailed calculation for the extended of SY inequality up to 1-loop and see the same 
expression with the original SY inequality applies in Sec. IV. In Sec V, we generalize the 
discussion in Sec. IV up to all order in the 5N expansion. In Sec VI, we summarize our 
conclusion. Detailed calculations in this paper are summarized in the Appendix A, B and C. 



2 Notations 



We explain some notations used in this paper. 

In order to express a tensor components, we use the lower-case alphabet indices such as 
a, 6, c . . . , running 1 to m. Then, we express a n-rank tensor as 



M aia2 .. Mn = M an = M n , 



(2.1) 



and the contraction for arbitrary two n-rank tensor M n and L n is expressed by using a metric 
G as , 



(M ni L n )gn =M an L an 



= g aibi g 



->a n b n 



a\a2- 



.a n -^bib2...b n 



(2.2) 



where (•, -)gn denotes an inner product made by the outer product of metrics, Q n = Q®Q®. . . , 
for n-rank tensors M n and L n . 

Furthermore, we define m(m n — l)/(m — l)-dimensional vectors M and L having n-rank 
tensors M n and L n as their components, 



L — (Li, L2, • • • , L n ) — (L ai Lab, • • • L aia2mmman ) , 

M = (Mi, M 2 , . . . , M n ) = (M 0J M abl M aia2 ... a J , 
and then we define the inner production for M and L as 

/go o\ ( M a \ 

M ab 



(2.3) 



(L, M)g = {L a , L a b, L a bc, 



J2(Ln,M n )gn 



g®g 
g®g®g 
\o .../ 



M abc 

V ••• / 



(2.4) 



where (•, -)g denotes a inner product defined by the metric G which is defined as 

g = g © {g ® g) e {g <g> g ® g®) e . . . 

'£0 

g®g 
.0 ... 



(2.5) 



The metric G is obviously positive-difinite when the metric g is positive-difinite. 
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3 The 5N formalism 



In this papar it is assumed that the universe is dominated by scalar fields in the inflationary 
period. In super horizon scales, the universe at each points behave as an independent uni- 
verse, called the separate universe. Furthermore, the primordial curvature perturbation £ is 
described by the fluctuation of the e-folding N = f* H'dt 1 ', where H is the Hubble parameter. 
Then, the curvature perturbation £ = SN can be expanded by the scalar fields at the initial 
time which is a few times of horizon crossing time t* as follows. 

C = 6N = N\p,tf(x)]-N\p,<fi] 

= ^—Sy a Jx) + - n h 6<p$(x)6<p b Jx) + . . . . (3.1) 
d(p% ; 2d(p%d(p\ ; ; v ; 

where the alphabets such as a, 6, . . . denote the number of scalar fields. Note that in the SN 
formalism, the e-folding must be expressed as a function of the background energy density p 
and scalar fields evaluated at initial time [9, 27-30]. Here, we express the SN expansion as, 

oo 1 

( = SN = J2- ] N a J<pi n , (3.2) 



71=1 



where we define as 



N = N 1 1 (3 41 



6<p? =^ 1 ^* 2 -"^* n , (3-3) 
d n N[p, gg i 

d^ 1 . . . dffi 

and N an is a complete symmetric tensor. 

In the calculation of the 5N formalism, we need the information about initial values of 
scalar fields. Therefore, we adopt a following assumption about initial conditions: 

• Fluctuations in the scalar fields at the horizon crossing S(p^ are scale invariant as well 
as Gaussian. 

The linear power spectrum for the scalar fields at initial time k (aH)* is assumed to be 

2tt 2 

(5^(Q6 V b Lk ,(Q) = ( 2 ^) 3 5 D (k + k0^3-ef , (3.5) 

where Stp^it*) denotes the Fourier transformation of scalar field perturbations. Since the 
linear perturbations of scalar fields have Gaussian distribution, the Non-Gaussianity of the 
curvature perturbation will arise from non-linear terms in the 5N expansion on super horizon 
scales. While usually the power spectrum Q® b is assumed to be proportional to 6 ab we do 
not adopt this assumption in this paper. This allow us to treat the case where two fields have 
a correlation each other at the initial time. The initial power spectrum Q® b is a real positive 
symmetric matrix and behaves as a positive metric in contracting for coefficients of the 8N 
formalism. In a diagrammatical picture of the 5N formalism, the non-dimensional initial 
power spectrum plays a role of the coupling constant and n-loop contributions contains 
the terms proportional to Q™. Then, the diagrams of the power spectrum, bispectrum and 
trispectrum are expressed by using the linear power spectrum as propagators (see Appendix 
A and B). 

Although the 5N formalism could be used for vector fields, we focus only on scalar 
fields in this paper, because the power spectrum for vector fields usually have strong scale- 
dependence and our proof in the below might not apply. The 5N formalism and the consis- 
tency relation for vector fields are studied in [34, 35]. 
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4 Tree level SY inequality 



As in the case of the original SY inequality, which is the most simple case, we first consider 
non-Gaissianity parameters Jnl and t^vl arising from the terms of the order 8ip 2 in the 8N 
expansion. Namely we consider the following form of the 5N expansion. 

C = N a 5vt + \N ab 5v a M. (4.1) 

This corresponds to the tree level calculation in a diagrammatical evaluation of the 6N for- 
malism. It has been shown that the following simple forms of the power spectrum, bispectrum 
and trispectrum are obtained in this case, respectively, 

V c = VfN a N b = N a N a , (4.2) 



iy a iy ab iy bc iy c , A A \ 

T NL = 7T T AM3 , (4.4) 



6 = N a N b N ab 
5 JNL (N C N C 

N a N ab N bc N c 

(N d N d y 

where the coefficients of the SN formalism, such as N a , N ab , etc., are contracted by using 
the power spectrum of scalar fields Q% b . 

Here, we use the Cauchy-Schwarz inequality. When we define the inner product for 
arbitrary two vectors v a and u a as (v,u) = Q^VaUb, then the Cauchy-Schwarz inequality leads 
to (v,u) 2 < (v,v)(u,u). Choosing v a and u a as v a = N a and u a = N^N^, we immediately 
find the original SY inequality. 

6 \ 2 = (N a ,N ab N b ) 2 
5 JNL J (N C N C ) 4 

. f (N a ,N a )(N bc N c ,N bd N d ) \ 

- { (w ) = TNL - (45) 



5 1-loop SY inequality 

In this section, we extend the SY inequality up to 1-loop order. This corresponds to the case 
where the 5N expansion is truncated at the order of 5(f 4 , namely, 

C =NJ<p a * + ^N ab 6<p a J^ 

+ j } N abc 8v a M h Ml + j,N abcd 5v a Mjv c M*. (5.1) 

Thus, we shall ignore the contributions in the power spectrum, bispectrum, or trispectrum 
coming from terms of 0(5(f 5 ). The 4th-order term is necessary when we calculate all of the 
1-loop contributions in /nl and tnl- Up to this order, not only all of the 1-loop contributions 
are included, but also some of the higher-order loop contributions are included. 
Then, the power spectrum is given, up to the 4th order, by 

V c (k) =N a N a + N ab N ab ln(feL) + N a N abb ln(/c max L) 

+ -^N acc N abb ln 2 (/c max L) + N abcc N ab ln(fcL) m(fc max L) . . . , (5.2) 
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where the 1st term is the tree contribution, the 2nd and 3rd terms are the 1-loop contribu- 
tions, and the 4th and 5th terms are the 2-loop contributions. The fc max is the ultra-violet 
cutoff. The L is a finite size of a box which is chosen to be much larger than the region of in- 
terest, namely L satisfies the condition Lk ^> 1 for any interest k. As pointed out in [36], the 
L is not large enough to satisfy ln(fcL) ^> 1, but it is the order of 1/Hq ~ ko = 0.002 Mpc -1 , 
where Hq is the current Hubble parameter and ko is the usual normalization scale used by 
the WMAP collaboration. 

This result can be simplified by defining the following quantities (see Eq. (25) of [37]): 

N a = N a + ^N abb ln(fe max L), (5.3) 
N ab = N ab + ^N abcc ln(fe max L), (5.4) 
N a bc = N abc + ^N abcdd ln(fc max L). (5.5) 



Then Eq. (5.2) becomes 



V c (k) = N a N a + N ab N ab In(feL) + . . . (5.6) 



Here, the dots in Eq. (5.6) include the higher-order terms such as N^ bcc . This is a nice way 
of writing the power spectrum etc., as the results do not include the ultra-violet cutoff, fc max , 
explicitly: the cutoff can be absorbed by redefining the derivatives of N. From now on, we 
shall remove the tildes from the equations, i.e., N — > N. 

As shown in Appendix A, the bispectrum and trispectrum with 1-loop corrections are 
calculated to be 

pN L (k bl , k b2 ) = V^\k bl )V^\k b2 ) [N a N ab N b + N a N abc N bc Qn(k bl L) + ln(fc &2 L)) + N ab N bc N ca In (k bl L) 

(5.7) 

TNL(k tl ,k t2 ) =V^\k tl )V^(k t2 ) [N a N ab N bc N c + 2N ab N abc N cd N d ln(k t2 L) 

+ {N ab N bc N cd N da + N a N abc N bcd N d + 2N a N abc N bd N dc ) ln(fc tl L)j , 

(5.8) 

where we have used pole approximation to evaluate the main contribution of the loop in- 
tegral(see Appendix A for more detail). The scale dependence is satisfied by the following 
conditions; k bl k b<2 and k tl k t2 . These expressions have been again simplified by using 
the redefinition of the derivatives of N and ignoring the higher-order terms. Here, we can 
set that k bl = k tl = k^ and k b2 — k t2 = ks without loss of generality, because the squeezed 
limit is satisfied for the bispectrum also in this setting. Therefore, we might focus only on 
the relation between fNL(kL,ks) and ttvl(^l? ks)- 

Now we are ready to prove the SY inequality. First we use the following inequality 
between arbitrary tensors such as M a: M ab and L ai L abl which is an application of the 
Cauchy-Schwarz inequality, 

(M a L a + M ab L ab f < (M a M a + M ab M ab ) (L a L a + L ab L ab ) , (5.9) 

where a positive-definite metric is used in contracting the indices of tensors. 
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Choosing M ai M ab and L ai L ab as follows. 

M a = N a , 

M ab = N ab [Hk L L)f\ 

L a = N ab N b + N abc N bc \n(k s L), 

Lab = (N abc N c + N ac N cb ) [\n{k L L)f 2 , 

Then it is easy to show the following inequality. 

'6 



(5.10) 



:fNL (k L ,k s ) 



< T NL (k L ,k s ) + 



N ab N abc N cde N de [Hk s L)Y 

r ( (k L )v 2 c (k s ) 



(5.11) 



where the last term in the right-hand-side (RHS) of Eq. (5.11) is the "2 loop" term. This 
result shows that, when we allow ourselves for completely general models in which this 
particular 2-loop term can become important, the SY inequality, tnl > (|/nl) > niay be 
violated badly. The reason of this violation is that we trancated the 5N expansion at a finite 
order. 

Still, from a model-building point of view, it is reasonable to assume that the 2-loop 
terms are sub-dominant compared to the tree or 1-loop terms; otherwise, we would have 
to require fine-tunings between the derivatives of N. Let us then study the consequence of 
ignoring this particular 2-loop term. We shall impose the following condition: 



N ah N ahc N cde N de [\ii{k s L)f 
N b N ba N ac N c 



N d N da N abc N bc ln(k s L) 



N b N ba N ac N c 



(5.12) 



This condition is ( 1-loop) » (2-loop) for the second term in the RHS of Eq. (5.8) and the 
last term in the RHS of Eq. (5.11). Interestingly, from the Cauchy-Schwarz inequality for 

v a N ab N b and u a = N abc N bc , we find 



/ N d N da N abc N bc \n(k s L) 
\ N b N ba N ac N c 

and from Eq. (5.12) we find 

N d N da N abc N bc ln(k s L) 



< 



N ab N abc N cde N de [\u{k s L)f 
N b N ba N ac N c 



(5.13) 



N b N ba N ac N c 



N d N da N abc N bc ln(k s L) 



N b N ba N ac N c 



(5.14) 



Finally, we obtain the following restriction on a particular form of 1-loop contributions: 

N d N da N abc N bc ln(k s L 



N b N ba N ac N c 



< 1. 



(5.15) 



As a result, if we ignore the last term in the RHS of Eq. (5.11), we must also ignore the 
second term in the RHS of Eq. (5.8). This is a peculiar feature of these terms, whose 
physical meaning is not clear. 

In any case, provided that the following additional condition is met: 
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• The 2-loop contributions are sub-dominant compared to the 1-loop contributions (or 
at least the particular 2-loop term in the RHS of Eq. (5.11) is small compared to the 
others), 

we finally arrive at the extended version of the SY inequality up to 1-loop: 

6 



rNh(k L ,k s ) > y-hh(k L ,k s )j , (5.16) 

which is valid as long as the 2-loop contributions are small. As pointed out in [38], we do 
not have the equality in the cases of single field models (single degree of freedom models) 
any more. 

It is thus proved that the same expression with the original SY inequality is satisfied 
even with 1-loop corrections. This result is a stronger than one in our previous work [32], 
and the previous conclusion that the 1-loop corrections weaken the original inequality should 
be modified. 

6 Generalization up to any loop SY inequality 

Now we prove the most general case where all order of the loop corrections are taken into 
account. The detail of the calculation for the power spectrum, bispectrum and trispectrum 
with any loop corrections will be given in the Appendix B. 
The power spectrum including all loop corrections is 

oo 

71=1 

OO 

= E© (Nn^gnQnkLr- 1 

71=1 

= £)(P„(fe), P n (k))gn = {P(k),P(k)) G , (6.1) 

77 

where we defined n-rank tensor P n depending on k as 

P n (k) = x [^(\nkL)^N n . (6.2) 
V n - 

The fjsiL and tnl with all loop contributions are given by 

~ oo oo oo 

SV 17 SV Z/ mi =1 7772=0 7773 = 1 

mims 



7Tli 17712 !77l3! 



oo oo oo oo oo oo 



_j W W 

TN L (k tl ,k t2 )= E E EEEH^r^-'H^r 1 ^- 2 

/ (IMJ/^MjJ mi =l m2 =l m3 =l m4 =0/ 1 =0i 2 =0 



mim 2 m 3 

ml , m2 , m3'm4 , ^l'^2 , am i-M™4 J ""l 6 '! 6 ^ JV &m 2 /i 2 c m3 ^c m3 e h d m4 , 

(6.4) 
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where the squeezed limit k^ = k\ <C k<i ~ ks = fc& 2 and the collinear limit k tl = |ki + 
ki2 1 <C k\ ~ ^2 — ks — k4 = kt 2 are taken, respectively. Even in this case, we can assume 
kbi — kti = and k\) 2 — k t2 = ks, where k^ <C ks, without loss of generality. 

Now, we are ready to derive the consistency relation between /atl(^l 5 ks) and t/vl(^L5 ks). 

First, we express f^L in terms of (•, For this purpose, we express //yl as 

~ 1 oo oo a— 1 

j/KLfe.fe) - n . fe)P((M E E E [m^r> w- 1 

s sv a=l mi=l m2=0 



mi (a — 7712) 
raj.! (a — m2)!m2! 



1 



OOOO a 



V c (k L )P c (k s 



E E EM^'M^r 1 " 1 



a=l mi=l m2=0 

mi (a — m,2) 



_ 1 / _ \ 1 ^ I Ami Omo Omn 

m\\(a — m2j!m2! 1 2 2 



T^a — mo c i 



2 c a — rw^^my 



00 00 en 



777,1777,2 



mi!(Q! — 777,2)!777,2! 



Oj-m-y ^a — rh-2 ^a-m2 £777,2 



c rn,2 ' 



(6.5) 



where a is defined as a = 777,3 + ^2- In the second equality, the upper limit of 777,2 is replaced 
from a — 1 into a, because for 777,2 = & the right-hand-side (RHS) of Eq. (6.5) becomes zero. 
We have used 777,2 = ol — 777,2 in the third equality. Then, we find 



00 00 a 



-fNL(k Ll k s ) 
5 



V ( (k L )P c (k s 



EE EP nfciL ] a_1 P n W mi - 1 



a=l mi=l m2=0 
777,1 — 777,2 



+ 



777,1777,2 



mi!(o! — m2)!m2! mi!(a — m2)\m2 



m|"m2 "rri2 L a-m2 



.Nr. 



00 00 a 



— - -— T V Y [lnkLLr-^lnksL]™ 1 - 1 

sv S «=1 7771 = 1 7772=0 

1 mitt 

2 mil (a — 777,2)! 777,2! 122 2 21 



(6.6) 



Choosing M a as 

M a (k Li k s ) = E E o m ,7' Z ,Vr» ^ L L)^(\nk s L)^- l N bm2am N amiCa _ m2 , 

Z — ' Z — ' Z 777,1 ! (a — 777,9 ! mo! 21 1 2 



7771 = 1 7772=0 



(6.7) 



and we find 



6 \ 1 °° 

^fNL(k L ,k s )J = Vc{kL)Vc{ks) E^fe)' M »^' fcg ))g Q 

(P(k L ),M(k L ,k s )) G , (6. 



V c (k L )V c (k s ) 
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where M a is a symmetrized completely because it is contracted by the complete symmetric 
tensor P a . 

Next, we can use the Cauchy-Shwarz inequality for the inner product (•, -)q because the 
inner product is positive-difinite, and we find 



^ ^271 !^2^ \ (P(k L ),P(kL)) G (M(k L , k s ),M(k L , k s )) G 



V{(k L )V*(k s ) 
1 

V c {k L )Vl{k s ) 

TNL(k L ,k s ) (6.9) 



(M(k L ,k s ),M(k L ,k s )) G 



where the last equality will be shown in Appendix C. 
We finally arrive at main result: 

r N L(k L ,ks)> (pNL{k L ,k s )\ • (6.10) 

This inequality has the same expression with the original SY inequality with only one differ- 
ence. Namely, the equality is no longer satisfied even in the case of single field models. 

Note that this inequality contains some errors. First, we use the pole approximation in 
the evaluation of the loop integrals. Exactly speaking, the power spectrum, bispectrum and 
trispectrum calculated by using the pole approximation is different from the truly evaluated 
quantities. In fact, the investigation of the validity of the pole approximation is an interesting 
problem. In [39], it is shown that the pole approximation is indeed valid for bispectrum with 
1-loop corrections in the leading order. Furthermore, there is also a problem about the 
infrared divergences which are eliminated by introducing the size L of the observed patch of 
the late universe [40, 41]. Second, theoretically predicted bispectrum and trispectrum with 
all loop corrections do not coincide with the local form bispectrum and trispectrum used in 
the observation. Therefore, we take the squeezed limit k^ = k\ <C /c2 — k% = ks for fjsfL and 
the collinear limit k^ = |ki + k^l <C k\ ~ &2 — ^3 — &4 = ks for tnl to make the calculated 
bispectrum and trispectrum into the local forms approximately. While the obtained limits 
of the bispectrum and trispectrum are approximately local form and have the maximum 
values at these limits as in the case of the local form bispectrum and trispectrum, we cannot 
completely remove some errors. 

Nevertheless, the inequality Eq. (6.10) would be still useful for focusing on the obser- 
vation of a complete violation of the inequality, because the errors are considered not to be 
large enough to change our inequality, i.e., 

TNL < (j>fNL) • (6-11) 

7 Conclusion 

We calculated the power spectrum, birpsctrum and trispectrum including the contributions 
of all higher order loop corrections for general multi-field models of inflationary scenario 
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using SN formalism. Then, we derived the following consistency relation between non-linear 
parameters under the assumption that the pole approximation dominates the higher order 
loop integral. 



This is exactly the same expression with the original SY inequality. However, there is one 
difference. While the equality is satisfied for single field models in the original SY inequality, 
we do not have the equality any more in the extended SY inequality. 

There are some errors for the extended SY inequality. The uncertainty causes the 
difference between theoretically calculated power spectrum, bispectrum, and trispectrum 
and the corresponding spectrums in the observation. We have to use the pole approximation 
to treat the loop integral analytically, and further take the squeezed limit and collinear limit 
to make the calculated bispectrum and trispectrum by using the pole approximation into the 
observed non-linear parameters f^L and t^l- Although the dominant contributions to the 
non-linear parameters would satisfy the SY inequality, these errors might slightly violate the 
SY inequality. 

Therefore, the most interesting case is the observation of a complete violation of the 
inequality, that is t^l *C (f/jvx) • Then the observation implies that a large class of inflation 
models are rejected for the mechanism to generate the observed fluctuations, provided that 
(1) scalar fields generate the primordial curvature perturbation; (2) the 6N formalism can 
be applied; (3) fluctuations in the scalar fields at the horizon crossing are scale invariant as 
well as Gaussian. 
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A Loop integrals: 1-loop 

In this appendix, we show the detail calculations of the 1-loop integrals. Here, the linear 
power spectrum of scalar fields, which is behave as a propagator, is given by 



where Q% is the non-dimensional power spectrum and behaves as a positive metric in con- 
tracting for coefficients of the 5N formalism. The diagrams for non-linear parameters f^L 
and tjsil is shown in Fig. (1) and Fig. (2), respectively. 

A.l power spectrum 

The power spectrum including a loop contribution is described as 




(7.1) 



(A.l) 
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Figure 1. Diagrams of the non-linear parameter f^L with 1-loop corrections 



where the subscript L on the integral denotes that the integrand is set equal to zero at 
p < L _1 . This setting is necessary in order to cut out a sphere around each of the singularities 
of the integrand. These singularities yield from the flatness of the linear power spectrum of 
scalar fields at initial time P*, diverging logarithmically at the singularities. This loop integral 
is analytically calculable, 
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k+i/L P* P^-k 2 ^ 
ln(fcL), 



(A.3) 



where we have used kL ^> 1. The UR divergence do not appear in this calculation. If we 
assume as ln(fcL) ~ 0(1) with L ~ fco, the last term in the second line in the RHS of Eq.(A.3) 
might not be neglected strictly, but we contain the term into the ln(fcL) as the uncertainty of 
the order unity. We can regard this result as the factor ln(fcL) yield from each singularities, 
such as p = and p = k. 
Finally, we can find 

V c (k) = N a N a + N ab N ab \n{kL). (A A) 

A. 2 bispectrum 

The bispectrum with 1-loop integrals is 
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(A.5) 



The second term in the RHS of Eq.(A.5) is not calculable analytically, but we assume that 
the main contributions in the integral yield from the singularities of the integrand as like the 
case of power spectrum, called the pole approximation. Therefore, this integral is calculated 
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as follows, 
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(A.6) 



where the upper limit of integral is determined as the next singular point. 

On the other hand, the third term in the RHS of Eq.(A.5) is calculable analytically as 
like the power spectrum, and we find 
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(A.7) 



Finally, we can find the bispectrum contained the 1-loop contributions as follows, 

7^/7 7 , \ _ N a N ab N h + NgN^N^ (ln(fciL) + m(fc 2 L)) + N ab N bc N ca In (Min{ki, k 2 }L) 
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(Ai 



This form does not coincide with the local form bispectrum, but this form has a maximum 
value in squeezed limit, hence we can correspond the 1-loop bispectrum to the local form 
Inl by imposing the squeezed limit, k bl = k\ <C fc 2 ~ &3 = fc& 2 . Then we find 

(3 

^fNL(k bl ,h 2 ) =V^ l {k bl )V^{k b2 ) [N a N ab N b + AU7V 6c 7V ca ln (fc 6l L) + N a N abc N bc (\n(k bl L) + ln(fc &2 L)) 

(A.9) 



The diagrams for this /jvx are shown in Fig. (1). 
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Figure 2. Diagrams of no n- linear parameter tjvl with 1-loop corrections. 
A. 3 trispectrum 

The trispectrum including 1-loop integrals is 
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(A.10) 



The second term in the RHS of Eq.(A.lO) is not calculable analytically as like the case of 
bispectrum, but this integral is evaluated by using the pole approximation as follows, 
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The third and forth terms are calculable analytically, and we find respectively, 
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(A.13) 



The fifth term is calculated by using the pole approximation, but we need a notice. In 
the last term of the RHS in (A. 10), if we set the pole around k4, we do not get the trispectrum 
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corresponding to t^l- When we extract only terms corresponding to t^l, we find 
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(A.14) 



fc 3 |ki+k 2 | 3 fc| 

Finally, we find the trispectrum including 1-loop contributions as follows, 

T^MMM) = V^ik^Wki + k 2 |)p ( r 1 (fc 4 ) 

N a N ab N bc N c + N ah N hc N cd N da In (Min{ki, k 4 , |k 4 + k 2 |) + N a N abc N bcd N d In (|k 4 + k 2 |L) 
+ N a N abc N bd N dc (In (Min{k 4 , |k 4 + k 2 |}L) + In (Min{k 1; |k 4 + k 2 |}L)) 



+ N ab N abc N cd N d (ln(hL) + \n(k 4 L)) P c (fcx)P c (|ki + k 2 |)P c (fe 4 ) 
+ 11 perm.. 

In order to identify the local form tnl, we impose the counter collinear limit, 
k tl = |ki + k 2 | < k t2 = k\ ~ k 2 ~ h 3 ~ fc 4 , 

and we find, 

T N L(k tl ,k t2 ) =p- 1 (fc tl )p- 2 (fc t2 ) [N a N ab N bc N c + 2N ab N abc N cd N d ln(k t2 L) 

+ (NabN^NcdNda + N a N abc N bcd N d + 2N a N ahc N bd N dc ) \n(k tl L) 



(A.15) 



(A.16) 



(A.17) 



where the diagrams for t^l with 1-loop corrections are shown in Fig. (2). 

These expressions for /tvl and t^vl m this Appendix are different from ones in [42] and 

[31]. 

B Loop integral : all loop 

We give the power spectrum, bispectrum and trispectrum with all loop corrections from a 
diagrammatical method. Each diagrams have coefficients and ln(fcL) factors. Therefore, the 
all loop expressions are given by 

(Coefficients ) x (Diagrams) x (ln(fcL) factors) , (B.l) 

number of diagrams 

where the coefficients are given by 

/ coefficients of \ / combinations of \ / coefficients arising from \ , . 
I the 6N expansion J I the lines of propagators J \ loop integrals. J ' 
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Figure 3. The power spectrum including n — 1 loops. The RHS of this diagram denotes briefly that 
the n lines propagate momentum k. 

B.l Power spectrum 

The two-point function for ( with all order is roughly expressed by 

00 00 TV TV 

(CO -EE -^T-^T <<V*" 1 ^*" 2 > • (B.3) 

77/1- 77/2- 

77,1 = 1 722 = 1 

Fig. (3) means that the two-point are connected by n lines of the propagator and these 
lines propagate the momentum k. We first find that the number of propagators must be 
ni — ^2 = ^5 and the number of the combination to connect each lines is n\. Therefore we 
find that Eq. (B.2) for power spectrum is 

/ 1 \ / 1 \ t IN / coefficients arising from\ , N 
UJUj X(n!)>< l loop integrals. )' ^ 

Next, the loop integral for the power spectrum is given by 

d 3 Pl . . . d 3 p n - ^ (B.5) 

|k-pi Vn\ 6 Pi Pt 

We can evaluate this integral by using the pole approximation. We assume that the main 
contribution of this integral come from around each poles such as pi = 0, P2 = 0, . . . , 
p n -i = and k = pi + • • • + p n _i, and they generate the factor ln(fcL). This means that 
only one line chosen in n lines propagates the momentum k, and the other n — 1 lines are 
integrated, generating the factor ln(fcL). 
Therefore the coefficient becomes 

-\ (-\ x (n!) x n = - (B.6) 
n\ \n\ n\ 



and the scale-dependence is given by [ln(/cL)] n 1 . 

Finally, we can derive the expression of the power spectrum with all loop corrections: 



00 
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where N an N an is defined as 

iVa„iVa„ = Q? bl G? b2 ■ ■ ■ QT bn N aia2 ... an N hlb2 ... bn . (Bi 
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mi 777-2 




Figure 4. Diagram for bispectrum including (mi + vfi2 + ms — 2) loops. 
B.2 Bispectrum 

The three-point function for £ with all order is roughly expressed by 

OO OO OO jyr jyr AT 

(ceo =EEE^f^f^f (5dr*s v z»5 V sr>). (b. 9 ) 

711 =1 77-2 = 1 723 = 1 

Fig. (4) means that each three points are connected by the mi, m2 and 7713 lines of 
propagators and the two lines for mi and 1712 propagate the momenta momentum ki and 
ki + k3 = — k2- For example, when mi — m2 — 1 and 7713 — 0, the diagram means the tree 
level one. First, we find that the mi, m2 and 7713 lines is generated from each points 774, 772 
and 773 respectively and the combinations to connect each lines are mi!, 7772! and 7773! as like 
the case of the power spectrum. Therefore, Eq.(B.2) for the bispectrum is 

/ 1 1 1 \ / nil 772! 773! \ / coefficients generated from\ 

\ni\ 77 2 ! 773! J \mi!m 2 ! mi!m 3 ! m 2 !m 3 ! J v m i- m 2- m 3-J ^ Y 00 p integrals. J 

1 / coefficients generated from 

mi!m2!m 3 ! \ loop integrals 

Next, we consider about the contributions of the loop integrals. Without loss of 
generality, we can consider that the all 7773 lines are integrated and generate the factor 
ln(Min{ki, k2}). On the other hand, the two lines in mi and 7772 propagate the momentum 
ki, — k2, and the mi — 1 and 7772 — 1 lines are integrated, generating the factor ln(fciL) and 
ln(&2£). Then, the coefficient generated from the loop integral is 77747712 . We apply the same 
discussion for the cases such as all mi and 7772 are integrated, generating the coefficients of 



(B.10) 
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m2m3 and mim3 respectively. According to above discussion, we find 
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To make the bispecturm correspond the local form non-linear parameter /nl, we take 
the squeezed limit k^ = k\ <C &2 — k% = k^ , and finally we find 
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where indeed the case of mi = 7773 = 1 and 7772 = means the tree level. 
B.3 Trispectrum 

In the case of the trispectrum, we need to consider Fig. (5). According to the similar discus- 
sion as like the case of the bispectrum, we can find that Eq. (B.2) for trispectrum becomes 



1 / coefficients generated from\ 

777i!7772!777 3 !7774!/i!/2! \ 1°°P integrals J ' 
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For the coefficients arising from loop integrals, there are twelve kinds of coefficients, that is, 
we choose three from mi, 7772, 7773, 7774, Zi, I2 except for Ziraira2, /177737774, /2^2^3 5 h^i^4 
which generate the other local form corresponding g^L- However, we eventually can express 
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2 7ft3 ~k 4 



Figure 5. Diagram of trispectrum corresponding to tnl with (mi + + ms + m4 + l\ + Z2 — 3) 
loops. 




Figure 6. Diagram of trispectrum corresponding to qnl with (mi + vrt2 + m3 + 777,4 + 'i + h — 3) 
loops. 



the trispectrum only in terms of vn\vn2vn^ as like the case of the bispectrum, and we find 
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When we take the collinear limit k tl = |ki + ki2 1 <C k t2 = k\ ~ &2 ~ &3 ~ &4, we finally 

find 
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(B.14) 

where the case of vn\ — rri2 — 1713 — 1 and 7714 = l\ = I2 = means the tree level. 
Similarly, we can derive the trispectrum corresponding to gj\fL from Fig. (6). 
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When we take the double squeezed limit k tl = fci ~ &2 <C A?3 — A?4 = k t2 , we can derive qnl 
with all loop corrections as follows: 

. cocococococo 
/b 7 m 1 =lm 2 =lm 3 =lm 4 =0Z 1 =0i 2 =0 

mim2m3 



mi ! m2 ! r«3 ! m4 ! 1 1 ! I2 ! ^ 2 ^ m4 am 1 ^ am i em 3 bm 2 ^ b m 2 fm 4 c h N Cl i em3 d, 2 

(B.15) 

C Some calculations 

In this appendix, we show the following equations in Eq. (6.9), 

r N L(k L ,ks) = 1 2 - (M(k L .ks),M(k L .ks)) G . (CI) 
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When we focus only on (M, M)g, we find 
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(C.2) 

where we define as M\ = Min{777,2, 777,4, a — 777,2, « ~ ^4}, M2 = Min{m2, 777,4} and M3 = 
Max{777,2, 777,4}. 

In order to consider Mi, M2 and M3 specifically, we decompose the range of 777,2 and 
777,4 as follows: (1) a — 2777,2; (2) < 777,2 < ol — 777,2; (3) ol — 777,2 < ^2 < a. The case (2) and 
(3) are equal when we replace 777,2 by 777,2 = ol — 777,2. The same discussion for 777,4 can be also 
applied. Therefore, the summations for 777,2 and 777,4 are decomposed as 
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where m is defined as m = 777,2 = ^4- 

The first terms of Eq. (C.3) means that M 1 = M 2 = M 3 = m = a/2. Then the RHS of 
(C.2) is 
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where we defined as 777, = m — /. This term is described in Fig. (5) as the case of 777,2 — 777,4 — fh 
and l\ — I2 — I- 
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The second term of Eq. (C.3) means that M\ — M2 — mi and M3 = 7B4 — a/2. Then 
the RHS of (C.2) is 
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where we defined as m 2 = m 2 — /, 7714 = 7714 — / and / 2 = 7774 — 77i 2 + I and l\ = Z. 

The third term of Eq. (C.3) means that M x = M 2 = M 3 = m. Then the RHS of (C.2) 
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7774 > 777 2 7771=1 7773 = 1 7712=0 /=0 

mim 3 (m 2 + m 4 + 2/) 
X mi!m 3 !m 2 !m 4 !(/!) 2 d ™4^ a ™i ^l^^a iV ^ 2 ^ c ™3 iVc ™3 e ^™ 4 

CO CO CO CO CO 

= E E E E E[ lnfciL ^ 2+ ' Tl4+2 '~ 1 [ lnfcsL ] mi+m3 ~ 1 

7774>7772 7711 = 1 7773 = 1 7772=0 /=0 

1 mim 3 (m 2 + m 4 + 2/) 
X 2 rai!ra 3 !m 2 !m 4 !(/!) 2 d ^4^ a ^i a ^i e ^^2 6 ^2^ c ^s c ^s e ^^ 4 

CO CO CO CO CO 

+ E E E E E[ ln ^ L ] 2 " + ' 1+ ' 2 " 1 [ ln ^ L r i+m3_1 

771=0 7771 = 1 7773 = 1 / 2 >U Zl=0 

l mim 3 (2m + /i + Z 2 ) v N N N (C ti\ 

X 2 mi!m 3 !(m!) 2 /i!/ 2 ! d ^^2 a ^i a ^i e ^i 6 ^ 6 ^^ 2 c ^3 iVc ^3 e ^i^' v u - b J 

where we defined as m 2 = m — I. 1714 = a — m — I and m = m 2 = 7714 and I — l\ — / 2 . The 
last equal in Eq. (C.6) is shown by the fact that the replacement of rn 2 <H> /1 and 7724 <H> / 2 
does not the form of the equation. In the explanation using Fig. (5), this implies that the 
both cases of 7774 > m 2 , Z 2 = h = I and 7714 — m 2 = 777, / 2 > Zi are equivalent. 

Finally, the forth term of Eq. (C.3) means that Mi — M 2 — m 2 and M 3 — 7714. Then 
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the RHS of (C.2) is 

OO CO OO CO CO 777,2 

E E E E E £ML£r 1 Ms£n i4ffl8 - 1 

« > 7714 7771=1 7773 = 1 7774 > 777 2 7772=0 /=0 

X 

mi!m3!(m2 — /)!/!(« — 7714 — Z) 1(7714 — ^2 + 0' 
x iV// f N h Nh f N a 

^cy. — m^ — ljl^rn-^ Q>m\ 67714 — m 2 + ^ ^2 — ^ ^2 — ^ ^ Cm 3 6777,4 — m 2 + ^ « — ^4 — ^ 

00 00 00 00 00 00 

E E E E E E[ lnfciL r 2+ ™ 4+ ' 2+ ' 1_1 [ lnfc s L r i+m3 ~ 2 

(7774 -7772 )>(/2-^l) m l = ! ^3 = 1 7772=0 l'2>h h=0 

2mim 3 (7n 2 + 7774 + Z 2 + ^l) v AT AT Af 

(C.7) 

where we defined as 7712 = m-2 — ^, "^4 = « — "74 — I, h = w-4 — ^2 + / arid li = Z. 
Adding Eq. (C.5) and Eq. (C.7), we find 

OO OO OO OO OO OO 

E E E E E E[ ln ^ L r 2+ ™ 4+/2+/1 + lnfc ^r i+m3 ~ 2 

7714 > 771 2 7771 = 1 7773 = 1 7772=0 h>h h=0 

m 1 m 3 (m 2 + 7774 + h + h) N N N N (r ^ 

X mi!m 3 !777 2 !m4!/i!/2! a ™i^2 d ™4 7Va ™i e 'i 6 ™ 2 ^^ ^-^J 

where this term means a diagram of 7774 > 7772 and I2 > h in Fig. (5). 

Adding Eq. (C.4) and Eq. (C.6) and Eq. (C.8) , we finally can arrive at 

(M(k L ,k s ),M(k L ,k s )) G 



V c (k L )V 2 c (k s ) 



^ 00 00 00 00 00 00 

vAkrWHk.) E E E EEEN^]™ 2 ^ 1 " 2 " 1 ^^]™ 1 ^" 2 

CV L) c, \ b ) mi= i m3= i m2 =0 7774=0 z 1= o z 2 =o 

1 mi 7773(7772 + 7774 + h + J 2 ) v AT AT AT 

X 4 mi!m2!m 3 !m 4 !/i!/2! iVa ^i^ 2 ^4 iVa ^i^i 6 ^2 iV ^2^2 c ^3 iVc ^3^i^4 



00 00 00 00 00 00 



= 1 2 E E E EE[ ln ^ L ] m2+m4+/1+/2_1 [ ln ^ L ] mi+m3_2 

V d k L)r c {ks) mi=1 ma=1 ma=1 m4=Q Zi=Q /2=Q 

m 1 m 2 m 3 

X mi!m2!m 3 !m 4 !/i!/2! arn ^ dm ± ^ a ™i e h b ™2 ^^h^m^ 1 ^ cm^e h dm A 
= r{k L ,k s ). (C.9) 
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